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Abstract 

Range trees are multidimensional binary trees which are used to perform d-dimensional or- 
thogonal range searching. In this technical report we study the implementation issues of range 
trees with fractional cascading, named layered range trees. We also document our implemen- 
tation of range trees with fractional cascading in C++ using STL and generic programming 
techniques. 

1 Introduction 

This project is an implementation of range trees with fractional cascading, named layered range 
trees in C++ using STL and generic programming techniques. Range trees are multidimensional 
binary trees which are used to perform d-dimensional orthogonal range searching. Range trees 
were discovered independently by several people including Bentley[lJ, who also discovered kd- 
trees and Lueker, who introduces the technique of fractional cascading for range trees [7j. An 
introduction to orthogonal range searching, range trees and fractional cascading can be found 
in [6l[9]. In [2] there is a presentation of a project of efficient implementations of range trees in 
2-3 dimensions including layered ones and some experimental results. 



2 Complexity issues 



The range trees answer a d-dimensional range query in time ©(log*^ n + k), where n is the whole 
set of points and k is the set of reported points. The construction time and the space the tree 
consume are O(nlog'^^^ n). Using fractional cascading we can be benefited by a logn factor in 
the last level of the tree and the resulting time complexity become O (log'^"^ n + k) . Intuitively, 
fractional cascading perform one binary search instead of two in the last level. The optimal 
solution to the orthogonal range search problem is due to Chazelle OH] who propose a structure 
with time complexity 0(log'^n + k) and 0(n(logn/loglogn]'^^^) space consumption, where c 
is a constant. 



3 Range trees in CGAL 

Although, CGAL library [3j provides some classes for range trees there is space for optimizations 
in that package [lOj. Firstly, there is a lack of recursive construction of d-dimensional range 
tree and the only way to construct a range tree of dimension d is to build a tree of dimension 1 
and then make this an associative range tree of a new one which will have dimension 2. Then 
one must build a tree of dimension 3 with this tree as an associative tree and this technique 
continues until the construction of the whole d-dimensional tree. 
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In addition to that, the package uses virtual functions, which increases the run time and 
finally there is no fractional cascading. 

The proposed approach uses nested templates for the representation of the d-dimensional 
range tree which is defined in compilation time. The dimension of the tree must be a con- 
stant and defined in the compilation time. In the last level a fractional cascading structure is 
constructed. 

For example a 4-dimensional range tree of size n with different kind of data at each layer is 
given by the following nested templated definition. 



Layered_range_tree <DataClass , 




Layered_ranj 


5e_tree <DataClass , 




Last_range_tree <DataClass> 




> 


> tree (n ) ; 





Note that for each layer I < d — 1 the same class Layered_range_tree is used. The last two 
layers, in which the fractional cascading is implemented, use the Last_range_tree class. The 
DataClass has the definitions of each layer's own data along with the comparison operators. 

4 Software implementation 

Essentially, the project was implemented using the C-|— |- language and the STL library [TT]. 
Concisely, the design uses methods from object oriented as well as the generic programming 
style. 

Representation. The trees are represented as STL vectors. The tree traversals are imple- 
mented using index arithmetic i.e. node's i parent is [V^J, the left, right child of 1 is 2i-|- 1 and 
2i + 2 respectively. This method is optimal for a full, static, binary tree and in our case the 
third is always hold. In order to have a full binary tree we replicate the last (biggest in the fist 
dimension) point and in the worst case we have a tree the half of which is useless with no effect 
to the time complexity (the replicated nodes would not be visited). In this project we are inter- 
ested in the static case of range trees but the design is sufficient for a dynamic implementation 
in which the tree nodes must also have some extra pointers. On the other hand, dynamization 
of the fractional cascading structure is not trivial [8] . 

Construction. For the construction of the tree we need to sort the input data with respect 
to the first coordinate and build recursively (top-down) the main tree in linear time. For the 
associative trees we don't have to sort the input data again. We build the associative trees in 
bottom-up manner. Every node merge the sorted lists of its children in linear time starting 
from the leaves which are trivially sorted. Note that this is essentially the same algorithm as 
merge-sort. 

Memory consumption. Even the asymptotic complexity of space stated above ensures 
that range tree needs a lot of memory. The only constraint in the number of dimensions of data 
is memory. Moreover, from the asymptotic complexity follows that with fixed memory there is 
a trade of between the number of data and number of dimensions. For example see 

I http://users.uoa.gr/~vfisikop/compgeom/Layered-Range-trees/examples/Layer ed-Range-tree-10.cpp| 

for an example of a range tree over 10-dimensional data. 
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5 Future work 



The most important feature for improvement is in the case that the tree is not fuh. The 
solution proposed in the previous section, seams very suboptimal and is not implemented yet. 
Another point is that the implementation don't handle yet the case that two points have equal 
coordinates. 
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